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ABSTRACT 



We present X-ray observations of two young energetic radio pulsars, PSRs B1046— 5^ 
and B1610— 50, and their surroundings, using archival data from the Advanced 
Satellite for Cosmology and Astrophysics (ASCA). 

The energetic pulsar PSR B1046— 58 is detected in X-rays with a significance 
of 4.5(T. The unabsorbed flux, estimated assuming a power-law spectrum and 
a neutral hydrogen column density Nh of 5 x 10^^ cm^^ is (2.5 ± 0.3) x 10^^'^ 
ergs cm~^ s~^ in the 2 — 10 keV band. Pulsed emission is not detected; the 
pulsed fraction is less than 31% at the 90% confidence level for a 50% duty 
cycle. We argue that the emission is best explained as originating from a pulsar- 
powered synchrotron nebula. The X-ray counterpart of the pulsar is the only 
hard source within the 95% error region of the previously unidentified 7-ray 
source 3EG J1048— 5840. This evidence supports the results of Kaspi et al., 
who in a companion paper, suggest that PSR B1046— 58 is the counterpart to 
3EG J1048-5840. 

X-ray emission from PSR B1610— 50 is not detected. Using similar assump- 
tions as above, the derived 30" upper limit for the unabsorbed 2 — 10 keV X-ray 
flux is 1.5 X 10^^^ ergs cm^^ s^^. We use the flux limit to estimate the pulsar's 
velocity to be less than ~170 km s^^, casting doubt on a previously reported 
association between PSR B1610— 50 and supernova remnant Kes 32. Kes 32 is 
detected, as is evident from the correlation between X-ray and radio emission. 
The ASCA images of PSR B1610— 50 are dominated by mirror-scattered emission 
from the X-ray-bright supernova remnant RCW 103, located 33' away. 

We find no evidence for extended emission around either pulsar, in contrast 
to previous reports of large nebulae surrounding both pulsars. 
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Subject headings: stars: neutron — pulsars: individual: (PSR B1046— 58, PSR 
B1610— 50) — supernova remnants — X-rays: stars — gamma-rays: star 
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1. Introduction 



X-ray observations of radio pulsars provide a powerful diagnostic of the energetics and 
emission mechanisms of rotation-powered neutron stars. As magnetic dipole braking slows 
the pulsar, it loses rotational kinetic energy at a rate E = An^IPP^^, where / is the moment 
of inertia of the neutron star, assumed to be 10^^ g cm^, and P is the rotation period. Though 
pulsars have traditionally been most easily studied at radio wavelengths, only a small fraction 
(10^^ to 10^^) of the "spin-down luminosity" E manifests itself as radio pulsations. Instead, 
it is believed that a significant fraction of the luminosity emerges as a relativistic wind of 
positrons and electrons. 

When this wind is confined by the surrounding medium, an observable synchrotron 
nebula results. Measurements of the morphology and spectrum are essential for determining 
the content and energy spectrum of the wind, probing the ambient density, and under- 
standing the shock acceleration mechanism. The growing list of observable X-ray emitting 
rotation-powered pulsars allows the study of the properties of the population as a whole. 
Young pulsars constitute a particularly interesting subset to investigate due to their large 
spin-down luminosities {^ 10^^ ergs s^^). 

PSR B1046— 58 was discovered by Johnston et al. (1992) during a 1400 MHz survey 
of the Galactic plane for radio pulsars. PSR B1046— 58 has E = 2.0 x 10^^ ergs s^^, char- 
acteristic age Tc = P/2P = 20 kyr, and properties reminiscent of the Vela pulsar. The 
dispersion measure (DM) toward the pulsar is 129 pc cm^^. Using the Taylor & Cordes 
(1993) DM-distance relationship, the distance to PSR B1046— 58 is estimated tohe d ^ 3.0 
kpc. The large E/And'^ value (a useful indicator of the detectability of its putative pulsar 
wind nebula) make it a strong candidate for X-ray emission. Kaspi et al. (1999) have found 
evidence for 7-ray pulsations from an EGRET point source spatially coincident with PSR 
B1046-58. 
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PSR B1610— 50 was also discovered in the Johnston et ah survey. PSR B1610— 50 is 
the fourth youngest Galactic pulsar {tc = 7.4 kyr) and has E = 1.6 x 10^^ ergs s~^ and 
DM = 586 pc cm"'^. PSR B1610— 50 occupies an interesting position in the age phase- 
space of known pulsars: its age is between those of the youngest pulsars {tc ^ 2 kyr), hke 
the Crab and PSR B0540— 69, and those in their "adolescence" (r^ ~ 10-20 kyr), like Vela 
and PSR B1046— 58. Caraveo (1993) proposed that PSR B1610— 50 is associated with the 
nearby supernova remnant Kes 32. However, using the DM-derived distance of 7.3 kpc and 
assuming the system's age is Tc, the pulsar's implied transverse velocity is vt > 2300((i/7.3) 
km s~^, significantly higher than the mean pulsar velocity (e.g. |Lyne fc Lorimer 1994| ). If it 
is traveling with so high a speed, one might expect strong ram-pressure confinement of its 
wind and a correspondingly high X-ray luminosity from the resultant synchrotron emission 
(see e.g. [Cheng 198^ ; [Wang, Li, fc Begelman 199^ ; Pinley, Srinivasan, fc Park 1996| ; [Romani 
ICordes, fc Yadigaroglu 1997| ; [Wang fc Gotthelf 199^) . 



Table 1 presents information for PSR B1046— 58 and PSR B1610— 50 including position, 
spin parameters, and derived quantities (i.e. spin-down luminosity and characteristic age). 

EDITOR: PLACE TABLE |T] HERE. 

Previous X-ray work on PSRs B1046— 58 and B1610— 50 has been limited in scope. 
Becker & Triimper (1997) published a ROSAT-hand (0.1-2.4 keV) luminosity for PSR B1046-5^ 
of 4 X 10^^ ergs s~^, but gave no further analysis. Kawai, Tamura, & Saito (1998) presented 
ASCA data on both pulsars, but restrict their attention to a single Gas Imaging Spectrom- 
eter (CIS) image of each pulsar. For both PSRs B1046— 58 and B1610— 50 they reported 
the detection of a large nebula (tens of arcminutes) associated with each pulsar. This pa- 
per undertakes a detailed analysis of the archival ASCA data, with an emphasis on image 
analysis. 
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2. Observations 

ASCA ( [I'anaka, Inoue, fc Holt 199^) observed PSR B1046-58 on 1994 January 27 and 
PSR B1610— 50 on 1994 March 25. We present an analysis of the data obtained from the 
public archive. For both observations, data were taken with all four imaging spectrometers, 
each in the focal plane of its own foil mirror: two Solid State Imaging Spectrometers (SIS- 
0, SIS-1) employing charge coupled devices (CCDs), and two Gas Imaging Spectrometers 
(GIS-2, GIS-3) employing gas scintillation proportional counters. These spectrometers offer 
moderate energy (~5%) and imaging (~2') resolution in their ~1 — 10 keV energy band-pass. 
The SIS has superior imaging and spectral capabilities, while the GIS has a higher effective 
area above ~2 keV and a greater net observation time than the SIS. To facilitate pulsation 
searches, GIS data were collected in the highest time resolution configuration (0.488 or 
3.906 ms depending on data acquisition rate). SIS data were acquired in 4 CCD mode with 
16 s integrations (PSR B1046-58) and 2 CCD mode with 8 s integrations (PSR B1610-50) 
using a combination of FAINT and BRIGHT modes (see the ASCA Data Analysis Guide for 
details). The data were filtered to exclude times of high background contamination using the 
standard "REV2" screening criteria. This rejects time intervals of South Atlantic Anomaly 
passages. Earth block, bright Earth limb in the field-of-view, and periods of high particle 
activity. The resulting effective observation times per single detector are 18 ks (GIS) and 
15 ks (SIS) for PSR B1046-58, and 11 ks (GIS) and 8.4 ks (SIS) for PSR B1610-50. 

3. PSR B1046-58 
3.1. Image Analysis 
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3.1.1. ASCA Data 

Flat-fielded images were generated by aligning and co-adding exposure-corrected images 
from the pairs of instruments. Exposure maps were generated with ASCAEXPO, ASCA 
software which uses the satellite aspect solution, instrument map (GIS), chip alignment, and 
hot pixel map (SIS) to determine the exposure time for each sky image pixel. The exposure 
correction was highly effective in removing the GIS instrumental structure due to the window 
support grid. Figure la displays the resultant smoothed broad-band (0.8 — 10 keV) image 
for the GIS. The image reveals emission confined to a slightly oval ~4' x 7' region, elongated 
along the direction parallel to declination. Though the statistics are limited, the emission 
appears to be concentrated in regions near the top and bottom of the oval. 

The pulsar location, determined by radio interferometric measurements made by Stap- 
pers et. al (1999), is marked by a cross and lies near the bottom of the emission region. The 
dashed square indicates the region shown in SIS images (Figures lb— Id). The two ellipses 
represent the 95% and 99% positional error boxes of the 7-ray source 3EG J1048— 5840 
(Hartman et al. 1999; see below for additional discussion). Figure lb displays the central 
region of the broad-band (0.4 — 10 keV) image for the SIS, while Figure Ic and Figure Id 
show the soft-band (0.4 — 2 keV) and hard-band (2 — 10 keV) images for the same region. 
A cross marks the location of PSR B1046— 58 and the arcs represent the errors ellipses for 
3EG J1048— 5840. Due to its superior spatial performance, the SIS resolves the smooth 
emission seen by the GIS into four point sources (hereafter called Src 1, Src 2, Src 3, and 
Src 4) possibly situated in a region of faint diffuse emission. 

The number of detected counts at the pulsar position is too small to fully resolve the 
familiar cross pattern of the X-ray telescope (XRT) PSF, as the morphology is dominated by 
Poisson fluctuations (see Hwang & Gotthelf 1997 §2.2 for a discussion of the significance of 
peaks in similarly processed images). To estimate the significance of the detections, we ignore 



the complexities of the ASCA point-spread function and compare the number of photons 
collected from a small aperture centered on the source with that from a 12' — 18' diameter 
concentric annulus. The relatively small number of counts available in the source region 
makes the approximation reasonable. Unfortunately, the close proximity of the four sources 
complicates this analysis. An optimally sizedQ 4' diameter aperture captures the majority of 
flux from the pulsar position as well as an undetermined amount of flux from neighboring 
sources, resulting in an overestimation of the pulsar's putative X-ray emission. Use of a 
smaller 2' diameter aperture eliminates the contamination from the neighbors but neglects 
the flux in the broad wings (~3') of the XRT PSF. Rather than artiflcially inflating the 
signiflcance of a detection, we employ the 2' diameter aperture with the understanding that 
our calculations may underestimate the strength of a source. Using the formalism outlined 
by Gotthelf & Kaspi (1998), we estimate the detection signiflcance (the signal-to-noise ratio, 
S/N) with an expression that accounts for both the source and background variance. 

Table 2 presents relevant information for the four sources detected by both SIS detectors, 
including source positions, count rates, detection signiflcance a, and hardness ratio H, where 
H = counts(2 - 10 keV)/counts(0.4 - 2 keV). Src 1 is relatively hard {H = 0.70 ± 0.30) and 
is strongly detected with a signiflcance of 4.5cr (44 background subtracted counts). Src 2 
is primarily soft {H = 0.55 ± 0.18) and is strongly detected with a signiflcance of 6. la (69 
background subtracted counts). Src 3 is hard {H = 1.55 ± 0.82) and has a signiflcance of 
3.9cr (38 background subtracted counts), while Src 4 is the hardest source {H = 2.15 ± 1.70) 
and has a signiflcance of 3.4o". In an attempt to reduce the influence of the diffuse emission 
and contamination from neighboring sources, centroid positions for Src 1, Src 3, and Src 4 



^Here, optimal refers to an aperture that maximizes the number of source counts captured 
within the extraction region, relative to the background contribution. See e.g. Gotthelf & 



Kaspi (1998). 
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were determined from analysis of the hard-band image in Figure Id and the position for 
Src 2 was determined from the soft-band image in Figure Ic. Combining the ~15" centroid 
uncertainty for each source with the ~20" revised pointing uncertainty ([Gotthelf 1996| ) for 
ASCA results in an overall source position uncertainty of order ~25". 

Src 1 lies 20" from the radio position of the pulsar, within the positional errors. As- 
suming approximately 4 sources per square degree with comparable flux to Src 1 (Gendreau, 
Barcons, & Fabian 1998) and a SIS spatial resolution of 3', we estimate the probability of a 
chance superposition of Src 1 with the pulsar's position to be of order 0.008. Src 2 lies 92" 
away from the radio position, making it extremely unlikely that it is the X-ray counterpart 
of the pulsar. 

Our images differ significantly from those produced by Kawai et al. (1998) from the 
same ASCA data. In particular, we do not find any evidence for a large (~10 — 20') nebula 
around the pulsar. Our reanalysis of the GIS data and analysis of the SIS data provide sup- 
port for emission from the pulsar (Src 1) that is unresolved by the ASCA PSF. Our disparate 
conclusions result primarily from the consideration of the SIS data with its superior spatial 
resolution, which shows that the oval shaped region in the GIS image actually represents 
emission from the four sources resolved by the SIS. We have also employed different proce- 
dures in the GIS analysis. Specifically, we have used an exposure correction that removes 
the significant structure produced by the GIS support grid and we have smoothed the data 
with a 3 X 3 boxcar function (~45" on a side) that approximates the core of the PSFQ 
When a much larger smoothing function is used, Poisson fluctuations, individual sources, 
and structure arising from the support grid (if not accounted for) can be blended into an 
apparent large, diffuse region of emission. 



^A similar procedure is also performed on the SIS data. The data is first rebinned x4, 
then smoothed with a 5 x 5 boxcar function (~30" on a side). 
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EDITOR: PLACE TABLE | HERE. 



3.1.2. RO SAT Data 

PSR B1046— 58 was also observed on 1996 March 8 by the HRI (High Resolution Imager) 
onboard ROSAT ioi 23 ks. Several point sources are clearly detected in the HRI FOV, but 
none is coincident with the position of the radio pulsar. The upper limit for a source at 
the radio position is <1 x 10"^ cps. This result is in disagreement with the detection of 
PSR B1046— 58 with the HRI reported by Becker & Triimper (1997) . A source is present, 
however, at 10M8'^13.''0, -58°30'44" (J2000), 80" away from the pulsar's radio position. 
This ROSAT source is only 19" away from SIS Src 2; given the apparently soft spectrum of 
Src 2 and the 25" positional uncertainty of the ASCA sources, suggests that they are the 
same source. 

To try to identify this source (SIS Src 2), we searched several optical catalogues. The 
only coincident source, located ~8" from the ROSAT position, was found in the Digital 
Sky Survey. Using the photometric calibration provided for the UK Schmidt Camera, we 
estimate a V magnitude of ~13.4. Positive identification of this optical source with the 
X-ray source detected by both ROSAT and ASCA requires spectroscopic data not presently 
available. 

There are no other ROSAT sources spatially coincident with the remaining three SIS 
sources, which is not surprising given their harder spectra. Furthermore, no extended or 
diffuse emission is seen in the HRI data. This suggests that the faint diffuse emission seen in 
the SIS data may not be physically significant. One possible origin of this diffuse structure 
could be a blending of Poisson fluctuations with emission from the closely grouped point 
sources scattered by the broad PSF of the ASCA mirrors. 
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3.2. Flux Estimation 

While the observations allow detection of the four sources, the low statistics prevent 
useful spectral analysis. However, the hardness ratio H of the SIS sources and the lack of 
i^OiS^T counterparts for all but Src 2 give some information about the sources. Srcs 1, 3, 
and 4 must either be absorbed non-thermal sources or thermal sources with temperature 
of at least several keV. Src 2 appears to be either an intrinsically soft thermal source with 
temperature on order of 50 eV or is a non-thermal source with a very steep power law that 
is undetectable above 2 keV in the SIS. 

To extract a flux estimate for the X-ray counterpart of the pulsar, we assume a spectral 
shape and adjust the overall normalization to match the count rate[| of Src 1. The canonical 
synchrotron nebula spectrum is characterized by a power law with photon index a = 2, 
where N{E) oc E^"" (see e.g. Reward fc Wang 1988| and Becker & Triimper 1997). The 



neutral hydrogen column density N^ can be constrained by combining the Seward & Wang 
(1988) approximation of 10 neutral hydrogen atoms per free electron with the DM or by 
using the HI maps of Dickey & Lockman (1990). The former yields Nh = 4.0 x 10^^ cm~^, 
the latter yields Nh < 1.4 x 10^^ cm~^. We adopt a value between these two rough estimates 
of Nh = 5 X 10^^ cm~^. After folding the spectral model through the instrument (SIS+XRT) 
response, we obtain an unabsorbed 2 — 10 keV flux of (2.5 ± 0.3) x 10^^'^ ergs cm~^ s~^. 
Folding these parameters through the i?05'74r instrument response yields an expected HRI 
count rate of (2.3 ± 0.3) x 10^'^ cps, in rough agreement with the upper limit (1 x 10^'^ cps) 
calculated above. The assumed spectral model and Nh also agree well with the observed 
hardness ratio. 



^Here, the count rate is that determined in §3.1, adjusted to account for the flux in the 
broad XRT wings that falls outside the extraction aperture. 
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3.3. Timing 

We carried out a timing analysis for PSR B1046— 58 using the combined data from the 
two GIS detectors. We selected events in the 2 — 10 keV band from a 4' diameter aperture 
centered on Src 1, using data acquired at the high and medium data rates only. A total 
of 472 events, a large fraction (~60%) of which are due to the background, were folded 
using an ephemeris obtained from radio timing observations of PSR B 1046— 58 at the 64-m 
Parkes radio telescope in New South Wales, Australia. Table 1 contains the ephemeris. As 
the putative pulse shape is unknown, we employed the H-test ( |de Jager 1994] ) to search for 
pulsations. For a duty cycle 6 = 0.5, the 3cr upper limit to the pulsed fraction is 0.31. For 
increasingly sharper pulse shapes oi 6 = 0.3 and 0.1, the 3a upper limits are 0.22 and 0.12, 
respectively. The absence of pulsations from the GIS data is consistent with the work of 
Saito (1998). 

4. PSR B1610-50 

4.1. Image Analysis 

Flat-fielded images were generated using the same prescription given in §3.1. Figure 2a 
displays the resultant smoothed broad-band (0.8 — 12 keV) image for the GIS. A cross marks 
the location of the pulsar determined from radio interferometric measurements (|Stappers 



Gaensler, fc Johnston 1999D n. The dashed rectangle shows the SIS FOV. The black contours 



are an overlay of 843 MHz MOST observations of the supernova remnant Kes 32 ( |Whiteoak 



fc Green 1996|) . The flux in the lower left quadrant results from scattered emission from 



^The interferometric pulsar position differs from the published catalog value by 57" (Tay- 
lor et al. 1995). No signal is present at either position. 
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the X-ray bright supernova remnant RCW 103, located 33' from the GIS optical axis. The 
scattered intensity has the gradient and shape expected when the 12' extent of the SNR 
is folded through the broad wings of the ASCA XRT QGotthelf, Petre, fc Hwang 1997| ). 
Examination of the soft-band (0.8 — 2 keV) and hard-band (2 — 12 keV) images in Figures 2b 
and 2c reveals that the contamination is largely confined to -E < 2 keV, due to the intrinsic 
spectral nature of RCW 103 and the decrease in the XRT scattering as a function of increasing 
energy. The SIS images have the same properties as the GIS images; for brevity, we only 
present the hard-band (2 — 10 keV) image in Figure 2d. 

We again ignore the complexities of the PSF and search for emission from the pulsar by 
comparing the number of photons collected from an (optimal) 4' diameter aperture centered 
on the radio location and with those collected from a 6' — 11' diameter concentric annulus. 
By restricting our search to E > 2 keV, we greatly decrease the amount of scattered emission 
from RCW 103. Our choice of background annulus avoids emission from Kes 32 and roughly 
contains the same amount of scattered flux as the source aperture region, allowing a reliable 
significance calculation. No emission was detected by either the GIS or SIS, and the combined 
detection significance is below 2a. 

As first noted by Kawai et al. (1998), the ASCA observation provides the first X-ray 
detection of the supernova remnant Kes 32. The low statistics prevent a detailed comparison 
of the X-ray emission with the elongated, shell-like radio morphology. To first order, the 
X-ray flux traces the radio intensity particularly along the western rim, as is evident in both 
the SIS and GIS images. The absence of X-ray emission from the direction of PSR B1610— 50 
contradicts a previous report of a large nebula powered by the pulsar ( |Kawai, 'lamura,~S 



Saito 1998| ). Contamination from the scattered RCW 103 emission and effects of smoothing 
with a function larger than the size of the ASCA PSF (see the discussion in § p.l.l| ) can 
account for the discrepancy. 
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Two additional sources are also present in the GIS data. At the top of the GIS FOV, the 
Einstein source 2E 1611.1—5018, a low-mass X-ray binary with J2000 coordinates 16'^14™54^, 
— 50°26'21", is clearly visible. This source is detected in both the hard and soft bands and 
was also detected by ROSAT. The second source, located south of the pulsar position and 
just outside the SIS FOV, is only seen above 2 keV. Data from both GIS detectors provide 
a 4.1(7 detection (47 background-subtracted counts). Its J2000 coordinates are 16'^14™18**, 
— 50°56'43", with a position uncertainty of ~1'. The source-like enhancements along the 
south-eastern edge of the GIS FOV result from scattered flux from RCW 103 and image 
processing artifacts. 

4.2. Flux Estimation 

The non-detection of PSR B1610— 50 can be used to place an upper limit on the flux 
from the pulsar. Starting with the observed background rates^ of 8.2 x 10~^ cps, we derive 
a 3a upper limit on the pulsar's count rate. We restrict our analysis to the GIS data, as 
these have a larger field of view for background estimation. We again use the canonical 
synchrotron nebula spectral model to estimate the flux and constrain the column density 
following the approach taken in §3.2. We adopt a value of Nh = 2 x 10^^ cm~^, consistent 
with the Seward & Wang estimate and the Dickey & Lockman upper limit. After folding 
the spectrum through the appropriate instrument response (GIS-I-XRT), we calculate a 3a 
upper limit to the unabsorbed 2 — 10 keV flux of 1.5 x 10~^^ ergs cm~^ s~^. 



^The background count rate is for an extraction region of radius 2'. 
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5. Discussion 

The importance of the detection of weak emission from PSR B1046— 58 and the non- 
detection of PSR B1610— 50 is most readily understandable in the context of the growing 
body of work on the X-ray properties of young (r^ = P/2P < 10^ yr) rotation-powered 
neutron stars. More than twenty of these objects have been detected, with three distinct 
physical processes responsible for the observed X-ray flux. 

Thermal emission can result either from the initial cooling of a young neutron star 
(e.g. [Ogelman 1995| and references therein) or from polar-cap reheating in older pulsars (e.g. 



Pecker fc Triimper 1993|) . If the heated cap sweeps across our line of sight, this emission may 
also be pulsed, as in the case of the Vela pulsar, which has an 11% pulsed fraction in the 
0.1—2.4 keV band ( Ogelman, I'lnley, fc Zimmermann 199iJ| ), or PSR B1055— 52, which has a 



similar pulsed fraction that varies with energy in the 0.1 — 10 keV band ([Ogelman &: Finley 



|1993| ; [Greiveldinger et al. 1996| ). The blackbody emission is independent of the spin-down 
luminosity E. 

Non-thermal magnetospheric emission, produced either in the polar cap or outer gap 
regions, is responsible for the classic "pulsar phenomenon," characterized by sharp pulsations 
of high pulsed fraction. The most famous example is the Crab pulsar, whose pulsed X- 
ray spectrum is characterized by a power law with photon index 2 ( P?oor fc Seward 1974| ). 
The energy for pulsed magnetospheric emission originates from the spindown and is seen 
from pulsars having high E at the extremes of the age distribution, from young [t^. < 
W yr) pulsars like the Crab, PSR B1509-58 and the two LMC pulsars PSR B0540-69 
and PSR J0537-6910, to old (r^ > 10^ yr) millisecond pulsars like PSR J0437-4715 and 
PSR B1821-24 ( ISeward fc Harnden 1982] ; [Seward, Harnden, fc Helfand 1984| ; |lVlarshall et 



al. 19981 ; [Becker fc 'Irlimper 1993| ; |Saito et al. 1997|) . 



The pulsar can have an associated synchrotron nebula, the observable evidence for the 
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channeling of the vast store of rotational kinetic energy into a wind of relativistic particles. 
Such nebulae arise when the pulsar wind is confined and subsequently shocked by the sur- 
rounding ISM, e.g. the nebulae produced by the Crab, Vela and PSR B1509— 58 ([Brinkmann 



Aschenbach, fc Langmeier 1985| and references therein; [Harnden et al. 1985| and Pgelnian fc 



Zimmermann 1989| ; Reward et al. 198^ , [Brazier fc Becker 19971 and references therein). 



5.1. PSR B1046-58 

For an age of 20 kyr, cooling models (see e.g. Ogelman 1995) predict thermal emission 
from PSR B1046— 58 to have an effective surface temperature of at most kT ^ 120 eV and 
a maximum bolometric luminosity of 2.3 x 10^^ ergs s~^. Assuming a 10 km neutron star 
radius, a 3 kpc distance and N^ ^ 5x 10^^ cm~^, the ASCA count rates should be no higher 
than 8 x 10^'^ cps (SIS) and 2 x 10^'^ cps (GIS). While these rates are comparable to the 
observed rates from the pulsar direction (refer to Table 2), the predicted count rate should 
fall to undetectable levels above 1.5 keV, in contradiction with the observations. Thus, 
cooling thermal emission cannot produce the observed flux from PSR B1046— 58. 

The apparently hard spectrum of the radiation suggests a non-thermal origin, either 
from the magnetosphere or from a synchrotron nebula. Magnetospheric emission is strongly 
pulsed, and given the upper limits on pulsations from the GIS data, it is extremely unlikely 
that magnetospheric emission contributes any significant fraction of the flux. Deeper obser- 
vations could reveal pulsations arising from either the magnetosphere or the modulation of 
thermal, surface emission. 

The most probable source of emission is synchrotron radiation powered by a relativistic 
pulsar wind. The most famous pulsar wind nebula surrounds the Crab pulsar ( [Rees fc Gunn 



T974t [Kennel fc Coroniti 198^ ; [Emmering fc Chevalier 1987[ ; [Gallant fc Arons 199^ . In 
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pulsar wind nebulae, the wind of relativistic electrons and positrons (and possibly heavy ions, 
Poshino et al. 1992| ) are confined, accelerated at the reverse shock, and radiate synchrotron 



emission. Using the fiux range calculated in §3.2 the 2 — 10 keV luminosity is Lx = (2.7 ± 
0.3) X 10^^ (i| g ergs s~^, where ds.o is the distance in units of 3.0 kpc. The conversion efficiency 
e of the spin-down luminosity E into ^S'C^-band emission is e = (1.3 ±0.1) x 10~^. 

The ~3' broad wings of the ASCA PSF severely limit detailed morphological studies 
of the synchrotron nebula. As the emission appears consistent with a point source, we 
can place a conservative 3' diameter upper limit on its angular extent. The limited spectral 
resolution also prevents us from considering whether the faint extended emission surrounding 
the four SIS sources could be the supernova remnant associated with PSR B1046— 58, a 
reasonable speculation given the pulsar's apparent youth. Ultimately, the nature of the 
extended emission and its possible relation to PSR B1046— 58 will only be resolved through 
observations with sufficiently high spatial and spectral resolution that would, for example, 
allow the detection of emission line features commonly found in other young supernova 
remnants. 

5.2. 3EG J1048-5840 

On a list ranked by E /d?, a parameter that has proven to be an excellent indicator of 7- 
ray detectability, PSR B1046— 58 is the ninth, with six of the eight sources higher being 7-ray 
pulsars. PSR B1046— 58 thus represents an excellent candidate for observable high-energy 
7-ray emission. Indeed, in a companion paper, Kaspi et al. (1999) suggest an association 
between PSR B1046— 58 and the unidentified high energy 7-ray source 3EG J1048— 5840 
based on the detection of significant pulsations from the 7-ray source at the radio pulsar 
period. 
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Identifying X-ray counterparts to unidentified EGRET sources is a useful way of sig- 
nificantly reducing the uncertainty in the position of the 7-ray source, on the assumption 
that the source spectrum extends into the X-ray band, true for both blazars and pulsars 
(e.g. [Kubo et al. 1998| ; Pecker fc Triimper 1997|) . In Figure la, we show the 95% and 99% 



confidence spatial contours of 3EG J1048— 5840 overlayed on the broad-band GIS image (R. 
Hartman, personal communication). Fortunately, the archival ASCA observation covers the 
entire 7-ray error box. The three detected ASCA sources discussed above (§3.1) are the only 
significant X-ray sources in the field; one or more of them is therefore probably the source 
of the 7-ray emission. 

The ASCA source nearest (within the 95% contour) the best-fit 3EG J1048-5840 po- 
sition is Src 1, which we have identified with PSR B1046— 58. It therefore represents the 
most likely counterpart to 3EG J1048— 5840, and supports the evidence presented by Kaspi 
et al. (1999). Src 2, outside the 95% contour but within the 99% region, was shown (§3.2) to 
have a soft spectrum, hence is unlikely to be the 7-ray source counterpart. The unidentified 
hard-spectrum Src 3 lies well outside the 95% contour but just within the 99% contour, so 
we cannot formally preclude its being the counterpart. However, were it the counterpart, 
it would most likely have to be a second young, energetic pulsar in the field, as it has no 
obvious radio counterpart, hence cannot be a blazar, since all known 7-ray emitting blazars 
are bright radio sources ( [Mattox et al. 199^ ). 



5.3. PSR B1610-50 

The non-detection of PSR B1610— 50 obviously prevents a study of the emission char- 
acteristics of the pulsar. However, derived upper limits are important for studying and 
understanding the X-ray properties of young rotation powered pulsars as a population, par- 
ticularly given this pulsar's place in age-space between the youngest Crab-like pulsars and 
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the older Vela-like pulsars. 

At an age r = 7.4 kyr, the Ogelman (1995) cooling model predicts thermal emission 
from the surface of the neutron star with a characteristic temperature kT ^ 130 eV. With 
the large column density towards the pulsar and its distance, though, the absorbed flux in 
the 2 — 10 keV band would be F^ = 2.1 x 10^^^ ergs s^^ cm^^, orders of magnitude below 
the ASCA detectability threshold. Pulsed magnetospheric emission is also expected to be 
present, given that pulsars with similar E exhibit pulsed non-thermal emission. However, as 
the total flux in the pulsed radiation is usually less than that from a synchrotron nebula, we 
only used the non-thermal X-ray emission expected from the nebula in our calculation of an 
upper flux limit in §3.3. For a flux limit F^^ < 1.5 x 10^^'^ ergs s~^ cm~^ and a distance of 
7.3 kpc, the 2 — 10 keV luminosity is L^ < 9.6 x lO'^^rf^g ergs s~^. The conversion efficiency 
e of the spin-down luminosity E into ASCA-hand emission is e < 6.1 x 10~^. 

If we assume that PSR B1610— 50 has an X-ray emitting synchrotron nebula, it is likely 
to be confined by ram pressure, as there is no evidence from X-ray or radio observations of 
a confining shell (i.e. SNR) around the pulsar ( [Green 1998| ). Relying on the work of Arons 
& Tavani (1993) and others, Gotthelf & Kaspi (1998 and references therein) show that the 
cooling efficiency for relativistic pairs of positrons and electrons is: 



e 



r = ^-«^'°"(aS5)'THL^)"^'looi;^"li5). (1) 



where tf and ts are the time scales for the pulsar wind fiow and synchrotron cooling, a 
is the ratio of magnetic energy flux to the kinetic energy flux of the wind (a ~ 0.005 for the 
Crab pulsar, [Kennel fc Coroniti 1984|) , p is the ambient density, v is the pulsar's velocity, and 
7 is the postshock pair Lorentz factor. ([Thompson et al. 1995[ ). The bolometric luminosity 



in synchrotron emission for the highest energy pairs is Lg ^ eE. If PSR B1610— 50 has wind 
properties similar to the Crab pulsar, the velocity must be <^ 170 km s^^ given the absence 
of X-ray emission. This velocity is inconsistent with an association between the pulsar and 
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Kes 32. Unless the pulsar resides in an extremely underdense region, i.e. p <^ 5.4 x 10"'^ 
cm~^, or has very different wind properties compared to the Crab pulsar, PSR B1610— 50 
does not have the transverse velocity suggested by Caraveo (1993). 

5.4. Lx — E relationships 

The luminosity measured for PSR B1046— 58 and the upper limit derived for PSR B1610— 50 
can be considered within the context of two previously determined empirical correlations 
between spin down luminosity E and X-ray luminosity L^: the Seward & Wang (1988) 
relationship (hereafter SW88) derived from Einstein data (0.2 — 4.0 keV band), logL^: = 
1.39 log -E — 16.6, and the Becker & Triimper (1997) relationship (hereafter BT97) derived 
from ROSAT data (0.1 — 2.4 keV band), L^. = lO^^E. To estimate the scatter associated 
with these models, we have computed the root mean square (RMS), defined here to be the 
square root of the variance, using the data presented by the authors in their papers. For the 
Einstein model, the RMS scatter is a factor of ~7, while for the ROSAT model, the RMS 
scatter is a factor of ~4. For PSR B1046— 58, SW88 predicts 7.2 ^-'q x 10^^ ergs s^"*^, compared 
to the observed (5.0±0.5) x 10^^ ergs s^^, while BT97 predicts 2.0q'5 x 10^^, compared to the 
observed (5.3 ±0.5) x 10'^^ ergs s^^. In both cases, the relationships overestimate L^ by fac- 
tors of at least several. For PSR B1610— 50, SW88 predicts 5.3q'^8 x 10'^^ ergs s^^, compared 
with the derived upper-limit 1.3 x 10^^ ergs s^^, while BT97 predicts 1.6 q'4 x 10^^, compared 
with the derived upper-limit 1.4 x 10^^ ergs s~^. For this pulsar, the derived upper-limits lie 
within the range of predicted values. 

One possible source of error is the distance derived from the Taylor & Cordes (1993) 
DM-d model. Distances are uncertain to ~25— 50%, translating to a potential error as a large 
as a factor of ~3 when converting flux to luminosity. Increasing the distance by ~25— 50% 
leads to a corresponding increase in L^ and helps improve the agreement between model and 
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observation. The discrepancies between model and observation also suggests that the current 
L,j.-E relationships, while certainly illustrating a correlation between spin-down and X-ray 
luminosity, may overlook important factors like the pulsar's velocity or the ISM conditions 
in its vicinity. 

6. Conclusions 

We have analyzed archival ASCA X-ray data for pulsars PSR B1046— 58 and PSR B1610— 50. 
We have detected emission from PSR B1046— 58 with a 4.5o" significance, the emission most 
likely arising from a synchrotron nebula. We estimate the 2 — 10 keV unabsorbed flux to be 
Fx = (2.5 ± 0.3) X 10^^^ ergs cm~^ s~^, assuming a distance d = 3.0 kpc. The data show no 
evidence for pulsation, with pulsed fraction less than 31% at the 90% confidence limit for 
a duty cycle of 50%. The synchrotron nebula is the only hard source within the 95% error 
ellipse of the previously unidentified 7-ray source 3EG J1048— 5840 and one of only two hard 
sources within the 99% error ellipse. The X-ray data, coupled with 7-ray pulsations discov- 
ered by Kaspi et al. (1999), strengthens the case for identifying PSR B1046— 58 as one of the 
few radio pulsars known to emit pulsed 7-rays. We detect no emission from PSR B1610— 50 
and place an upper limit on the 2 — 10 keV unabsorbed flux F^ < 1.5 x 10^^'^ ergs cm^^ s^^. 
This flux limit argues that the pulsar's velocity f ^ 170 km s^^, providing evidence against 
the association between by PSR B1610— 50 and Kes 32 claimed by Caraveo (1993). 

The flux-derived luminosities for both pulsars are below those predicted by the Se- 
ward & Wang (1988) and Becker & Triimper (1997) L^ — E relationship, particularly for 
PSR B1046— 58. The differences may represent the intrinsic limitations of the current em- 
pirically derived L^ — E relationships, which do not consider parameters such as ambient 
density and pulsar velocity. 
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The unresolved nebula powered by PSR B1046— 58 and the non-detection from PSR B1610— 50 
contradict previous reports by Kawai et al. (1998) of large (~10 — 20') nebula associated 
with these pulsars. These observations also demonstrate the difficulties associated with ana- 
lyzing and interpreting current hard X-ray {E > 2 keV) imaging data. The upcoming launch 
of the Chandra X-ray Observatory, with its unparalleled resolution (~0".5) and hard-band 
sensitivity, and XMM with its good spatial resolution (~8") and large collecting area, should 
improve this situation considerably. 

We thank B. Gaensler for providing the radio images of the PSR B1046— 58 and PSR B1610— 50- 
Kes 32 fields and for useful discussion. We thank M. Bailes, R. Manchester, and R. Pace 
for acquiring radio timing data for PSR B1046— 58 at the 64-m Parkes radio telescope in 
New South Wales, Australia. We thank R. Hartman for providing information on 3EG 
J 1048— 5840 prior to publication. We thank W. Becker for discussion on the analysis of 
the ROSAT observations of PSR B1046— 58. Finally, we thank the anonymous referee for 
useful comments and suggestions that improved this paper. We have used the NASA main- 
tained HEASARC web site extensively for archival data retrieval and subsequent analysis. 
M. Pivovaroff is supported in part by NASA under Contract NASA-37716. 



-23- 

REFERENCES 

Arons, J. & Tavani, M. 1993, ApJ, 403, 249 

Becker, W. & Triimper, J. 1993, Nature, 365, 528 

Becker, W. & Triimper, J. 1997, A&A, 326, 682 

Brazier, K. T. S. & Becker, W. 1997, MNRAS, 284, 335 

Brinkniann, W., Aschenbach, B., & Langmeier, A. 1985, Nature, 313, 662 

Caraveo, R A. 1993, ApJ, 415, LI 11 

Cheng, A. F. 1983, ApJ, 275, 790 

de Jager, O. C. 1994, ApJ, 436, 239 

Dickey, J. M. & Lockman, F. J. 1990, ARAA, 28, 215 

Emmering, R. T. & Chevalier, R. A. 1987, ApJ, 321, 334 

Finley, J. R, Srinivasan, R., & Park, S. 1996, ApJ, 466, 938 

Gallant, Y. A. & Arons, J. 1994, ApJ, 435, 230 

Gendreau, K. C, Barcons, X., & Fabian, A. C. 1998, MNRAS, 297, 41 

Gotthelf, E. V. 1996, Improved ASCA Sky Coordinates, (Greenbelt: NASA Goddard Space 
Flight Center), ( http://legacy.gsfc.nasa.gov/docs/asca/coord/updatecoord.html ) 



Gotthelf, E. V. & Kaspi, V. M. 1998, ApJ, 497, L29 
Gotthelf, E. v., Petre, R., & Hwang, U. 1997, ApJ, 487, LI 75 



-24- 

Green, D. A. 1998, A Catalogue of Galactic Supernova Remnants (1998 

September Version), (Cambridge: Mullard Radio Astronomy Observatory). 
( |http://www.mrao.cam.ac.uk/surveys/snrs/|) 

Greiveldinger, C. et al. 1996, ApJ, 465, L35 

Harnden, F. R., Grant, R D., Seward, F. D., & Kalin, S. M. 1985, ApJ, 299, 828 

Hoshino, M., Arons, J., Gallant, Y. A., & Langdon, A. B. 1992, ApJ, 390, 454 

Hwang, U. & Gotthelf, E. V. 1997, ApJ, 475, 665 

Johnston, S., Lyne, A. G., Manchester, R. N., Kniffen, D. A., D'Amico, N., Lim, J., & 
Ashworth, M. 1992, MNRAS, 255, 401 

Kaspi, V. M., Lackey, J., Mattox, J., Manchester, R. N., & Bailes, M. 1999, ApJ. in press 

Kawai, N., Tamura, K., & Saito, Y. 1998, Adv. Space Res., 21, 213 

Kennel, C. F. & Coroniti, F. V. 1984, ApJ, 283, 710 

Kubo, H., Takahashi, T., Madejski, G., Tashiro, M., Makino, F., Inoue, S., & Takahara, F. 
1998, ApJ, 504, 693 

Lyne, A. G. & Lorimer, D. R. 1994, Nature, 369, 127 

Marshall, F. E., Gotthelf, E. V., Zhang, W., Middleditch, J., & Wang, Q. D. 1998, ApJ, 
499, L179 

Mattox, J. R., Schachter, J., Molnar, L., Hartman, R. C, & Patnaik, A. R. 1997, ApJ, 481, 
95 

Ogelman, H. 1995, in The Lives of the Neutron Stars (NATO ASI Series), ed. A. Alpar, U 
Kiziloglu, & J. van Paradijs, (Dordrecht: Kluwer), 101 



-25- 

Ogelman, H. & Finley, J. P. 1993, ApJ, 413, L31 

Ogelman, H., Finley, J. P., & Zimmermann, H. U. 1993, Nature, 361, 136 

Ogelman, H. & Zimmermann, H.-U. 1989, A&A, 214, 179 

Rees, M. J. & Gunn, J. E. 1974, MNRAS, 167, 1 

Romani, R. W., Cordes, J. M., & Yadigaroglu, I.-A. 1997, ApJ, 484, L137 

Saito, Y. 1998. PhD thesis. University of Tokyo 

Saito, Y., Kawai, N., Kamae, T., Shibata, S., & Dotani, T. 1997, ApJ, 477, 37 

Seward, F. D. & Harnden, F. R. 1982, ApJ, 256, L45 

Seward, F. D., Harnden, F. R., & Helfand, D. J. 1984, ApJ, 287, L19 

Seward, F. D., Harnden, F. R., Szymkowiak, A., & Swank, J. 1984, ApJ, 281, 650 

Seward, F. D. & Wang, Z.-U. 1988, ApJ, 332, 199 

Stappers, B. W., Gaensler, B. M., & Johnston, S. 1999. in preparation 

Tanaka, Y., Inoue, H., & Holt, S. S. 1994, PASJ, 46, L37 

Taylor, J. H. & Cordes, J. M. 1993, ApJ, 411, 674 

Taylor, J. H., Manchester, R. N., Lyne, A. G., & Camilo, F. 1995, Unpublished (available 
at [ftp: //pulsar. princeton.edu/pub/catalogl) 



Thompson, D. J. et al. 1995, ApJS, 101, 259 
Toor, A. & Seward, F. D. 1974, AJ, 79, 995 
Wang, Q. D. & Gotthelf, E. V. 1998, ApJ, 494, 623 



-26- 

Wang, Q. D., Li, Z.-Y., & Begelman, M. C. 1993, Nature, 364, 127 
Whiteoak, J. B. Z. & Green, A. J. 1996, A&AS, 118, 329 



This manuscript was prepared with the AAS h^T[T;X macros v4.0. 



-27- 



Table 1. Astrometric and Spin Parameters for PSRs B1046— 58 and B1610— 50. 



Parameter 



PSR B1046-58 



PSR B1610-50 



Right ascension (J2000)^ 

Declination (J2000)^ 

Period, P (s)^'= 

Period derivative, P (s s~^)'^ 

Epoch of period (MJD) 

Dispersion measure, DM (pc cm~^) 

Characteristic age, Tc (kyr) 

Spin-down luminosity, E (ergs s^^) 

Distance, d (kpc)'^ 

Column density, A''^^ (cm~^)*^ 



10'"48°^12."6 



-58°32'03."8 

0.124 

95.9 X 10-15 

49403.0 

129.09(1) 

20.4 

2.0 X 10^^ 

3.0 

(0.40 - 1.4) X 10^2 



16hl4miisg 

-50°48'01."9 

0.232 

493 X 10-15 

48658.0 

586(5) 

7.4 

1.6 X 10^^ 

7.3 

ri.8-2.2) X 10^2 



'^Radio positions, from Stappers, Gaensler, Sz Johnston 1999 , are uncertain 

by < 0.1". 

'^Data obtained from the Taylor et al. (1995) pulsar catalog. This information 

is not used for the timing analysis of PSR B1046— 58. 

^Ephemeris for PSR B1046— 58 was obtained from radio timing observations 

at the 64-m Parkes radio telescope in New South Wales, Australia. 

^^Derived from the Taylor & Cordes (1993) DM-distance model. 

*^Lower limits were derived from the Seward & Wang (1988) estimate of 10 
neutral hydrogen atoms per free electron; upper limits were derived from Dickey 
& Lockman (1990). 
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Table 2. ASCA SIS detection of PSR B1046-58 (0.4-10 keV Band). 

Source Right ascension Declination Count Rate Background Rate Hardness Significance 
(J2000) (J2000) (xlO-^cps) (xlO'^ cps) Ratio, H (a) 



Src 1 


10'^48"^11.^6 


-58°31'46" 


3.99±0.42 


2.10 ±0.04 


0.70 ± 0.30 


4.5 


Src2 


10''48'"15.n 


-58°30'34" 


4.87±0.43 


2.18 ±0.04 


0.55 ±0.18 


6.1 


Src 3 


10''48'°04.M 


-58°28'26" 


3.62±0.38 


2.15 ±0.05 


1.48 ±0.82 


3.9 


Src 4 


10''48™15.n 


-58°28'57" 


3.27±0.38 


2.21 ±0.04 


2.15 ±1.70 


3.4^ 



Note. — The positions for Src 1, Src 3, and Src 4 were derived from the hard band (2 — 10 keV) image; 

the position for Src 2 was derived from the soft band (0.4 — 2 keV) image. Total uncertainties in the source 

positions are --^25". Count rate is the total source plus background count rate in an aperture centered on the 

source position. Background rate is the count rate in a 12' — 18' diameter annulus concentric with the source 

position, normalized to the source aperture. Refer to the text (§3 & §4) for the definition of significance 

and further discussion on all the measured quantities. The hardness ratio H is defined as counts in the hard 

band (2 — 10 keV) divided by counts in the soft band (0.4 — 2 keV). 

'^This significance is for the 2 — 10 keV band. The broad-band significance is 2>.2a. 
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Fig. 1 - ASCA images of the PSR B1046— 58 field: flat-fielded images of the region around the 
pulsar, whose location is marked by the cross, a) The broad band (0.8 — 12 keV) GIS image shows 
an oval shaped region of X-ray emission with the pulsar located at its southern tip. The two ellipses 
represent the 95% and 99% error boxes for the 7-ray source 3EG J 1048-5840. The dashed square 
delineates the SIS region displayed in b)— d). b) The broad band (0.4 — 10 keV) SIS image clearly 
showing the three labeled sources embedded in a diffuse emission region, c) The soft band SIS 
image (0.4 — 2 keV) revealing the soft, probably thermal nature of Src 2. Note that Srcs 1, 3, and 
4 are very weak in this band, d) The hard band (2 — 10 keV) SIS image showing the hard nature of 
Srcs 1, 3 and 4. We identify Src 1, offset 20" from the radio position of PSR B1046— 58 and the only 
source inside the 95% error circle of the pulsed 7-ray source 3EG J1048— 5840, as the synchrotron 
nebula of PSR B1046— 58. Contours approximately correspond to the 4o", 5a, 6a, 7a, 8a, and 9a 
levels. Count rates are in units of 10^^ cps/pixel for the GIS and 10~® cps/pixel for the SIS. 
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Fig. 2 - ASCA images of the PSR B1610— 50 field: flat-fielded images of the region around the 
pulsar, whose location is marked by the cross. The dashed rectangle represents the SIS FOV, 
and the contours, corresponding to 0.04, 0.18, 0.31, 0.44, 0.57, and 0.70 Jy bcam^^, are from 843 
MHz MOST observations of the supernova remnant Kes 32. a) The broad band (0.8 — 12 keV) 
GIS image of the PSR B1046— 58 field. Scattered emission from the nearby supernova remnant 
RCW 103 is responsible for the large flux gradient that begins in the southeast FOV and extends 
to th edge of the SIS FOV. b) The hard band (2 — 12 keV) GIS image shows enhanced emission 
that traces the radio emission from Kes 32. No significant fiux is seen from the pulsar location. 
The previously known Einstein source 2E 1611.1—5018 is visible at the top of the FOV, while an 
unidentified source is located approximately due south of the pulsar position, c) The soft band 
(0.8 — 2 keV) GIS image explicitly shows the extent of the scattered emission from RCW 103. Note 
the distinct lack of emission from Kes 32. d) The hard band (2 — 10 keV) SIS image similarly 
shows the correspondence between the radio contours and X-ray emission and no emission from 
PSR B1610— 50. Count rates are in units of 10"^ cps/pixel for the GIS and 10^^ cps/pixel for the 
SIS. 



